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Abstract 

When a hydrogen atom trapped in finite space formed by two parallel per- 
fectly conducting plates, the transition rates can be modified. Life-time of 
the hydrogen 2P 1( / 2 -state is estimated to be shorter as much as 56.3ps with a 
separation distance of b = 1.2//m between two plates. Although the modifi- 
cation is dependent on position of the atom, the life-time can be shorten or 
lengthen by adjusting the separation distance even for the case that the atom 
is placed at the center between two plates. 
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§1. INTRODUCTION 



A long time ago, Casimirfl] found that the zero-point energy appearing in quantization 
of the classical radiation field in free space could be observable by introducing two parallel 
perfectly conducting plates. He derived 



where L 2 is the area of plates and b is the separation distance of two plates. The first and 
second terms on the left hand side denote the zero-point energies of radiation fields quantized 
in finite and free spaces, respectively. Notice the non-vanishing value on the right-hand side 
which is inversely proportional to cube of b. This result implies spatial dependence of the 
zero-point energy. On the other hand, it was also shown that the atomic energy levels could 
be modified when the atom was placed between two parallel plates [2,3,4]. Generally, the 
energy levels are expressed by real parts of the energy eigen-values while the life-times of 
energy levels, equivalently transition rates, are given by imaginary values of the complex 
energy eigen-values. Therefore, one may expect that the atomic transition rates will also 
be modified as long as the energy levels are shifted in finite space[5,6]. In this paper, we 
investigate the atomic transition rate in the framework similar to Casimir's calculation of 
the attractive force exerting between two plates [1]. 



Let us start to briefly review the calculation of transition rate of the hydrogen atom in 
free space. 

The time-dependent perturbation theory in quantum mechanics gives the first-order am- 
plitude of spontaneous radiation as [7] 




(1) 



§2. CALCULATION OF THE TRANSITION RATE 



2.1 Case of free space 
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where the interaction Hamiltonian is expressed as 



H I (t') = H' I exv(±iut t ). (3) 
Then, the transition probability is 

\a^\t)\^ = ^\HH'M 2 \&, (4) 

where 

I$= [ t dt'exp[t(u ±u)t'] (5) 
Jo 

with Tiuq = E m — E x . The negative (positive) sign in eq.(5) denotes the photon emis- 
sion(absorption). Hereafter, we consider only the photon emission, i.e. negative sign 
in eq.(5). Carrying out the integration directly, we obtain 

( _) 2 = exp{i(uj - u)t} - 1 2 

= 4 Sin ;^°-^. (6) 
(u - id) 2 



Then, the transition rate is given by [7,8] 

h 2 J IN ' 11 /l (u; -u}) (2tt) 3 c 3 V ; 

If the dipole approximation is taken for Hj, eq.(7) can be given as 

ft 2 2\/J Y 1 71 (coo-co) (2ttc) 3 

e 2 /",, , ... |2 jri /°° ^ 3 sin[t(cu - u)] 



87T3RC3 7 l(m|r|/) • £| ^ loo (C -C) ^ 

where e is the polarization vector and d = er is a dipole moment. Considering that 

(m|r|Z) • e = r mi cos©, (8) 

and [7] 
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J cos 2 Qdn = y , (9) 

we find 

e 2 2 [°° w 3 sin[t(w - w)] , 



3lT 2 ?lC 3 J-oo (uJo — 

(/-) \4w ' 

\Anfic J 3 c 2 



r mi | 2 - (10) 



The matrix element r m \ can be calculated with the hydrogen wave functions and the result 
for the 1P\ji — > 15*1/2 transition is known as 

where ao is the Bohr radius. Since v{2S\/2 — > lSi/2) = 2466061395. 6(4. 8)Mifz[9] and the 
Lamb shift yields v(2S 1/2 -> 2P 1/2 ) = 1057.8594 ± 0.0019M#z [10], we find z/(2P 1/2 -> 
151/2) = ^o(2Pi/ 2 -> 15i /2 )/27r = 2466060337. 7Mi7z and, then, the transition rate is 

W / "o(2P 1/2 -> l.Si/2) = 0.626288 x 10 9 ifz, (12) 

which gives the life-time of the 2Pi/ 2 -state 

r (2Pi /2 ) = [w^o(2Pi/2 -> 15i /2 )] _1 = 1.597 x 10~ 9 s. (13) 

This is the result for the hydrogen atom in free space. 



2.2 Coordinate system of dipole moment and polarization vector 

Although the polarization vector, e, is perpendicular to the wave vector, k, of the radi- 
ation field, i.e. e • k — 0, the dipole moment, d, is generally not specified in any direction. 
It is better to clarify relations among these vectors. 

First of all, the space-fixed coordinates are assigned by unit vectors (ei,e 2 , e 3 ) and the 
wave vector of radiation field, k, is chosen in the direction of a unit vector e 3 . Here, unit 
vectors (e 1: e 2 , e 3 ) are orthogonal to each other. Then, we have 
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e 3 = e.\ sin 9 cos + e 2 sin # sin + e3 cos 9 (15a) 
where # is the angle vetween e 3 and e 3 and is the angle between e.\ and a projection of 



e 3 on the ei — e 2 plane. Since e x is perpendicular to e 3 , replacing # by # + | in eq.( |15a| ), 
we find 

e.\ = ei cos 9 cos + e 2 cos#sin0 — e 3 sin#. (15b) 

Considering that e 2 is perpendicular to the ei — e 3 plane, we replace by + | and set 
# = | in eq.( |15a| ). Then, we can have 



e 2 = — ei sin0 + e 2 cos0. (15c) 

Polarization of the radiation field is perpendicular to k, i.e. , e ■ k = and, then, it is 
expressed as 

* = e { (i= 1, 2). (16) 



If the dipole moment is along e 3 , i.e. d = de 3 , we find by eqs.( |15b| ) and ( |15c 



p2tt p2tt ( 

J | < m\{e ■ d)\l > \ 2 d<f) = J d 2 ml sin 2 9d<P = 2*^(1 - cos 2 9) = 2nd 2 ml \l-fA. (17) 
This result yields that given in eq.(4.9) of ref.ll. On the other hand, for d = de\, 

/*27T /•2-7T 

/ | < m|(e ■ d)\l > \ 2 d(j) = d 2 ml / [(cos 9 cos 0) 2 + (sin 0) 2 ]<i0 

JO JO 



2' 



= ttgQI + cos 2 9) = nd 2 ml (l + ^J, (18) 

which yields the result given in eq.(4.17) of ref.ll. For d = de 2 , we obtain the same result 
as eq.(18), 

j Q | < m|(e • d)|Z > | 2 rf0 = 7rd^(l + cos 2 9) = nd 2 rd f 1 + J . (19) 

Let us explore a general case that d is not polarized along a specific direction, i.e. 
d{9' , 0'). When the normal mode of the electric field is expressed Uj (j = x, y, z), we have 
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I < m\(e ■ d)u\l > | 2 = | < m\e ix d x u x + eiyd y d y + e iz d z u z \l > \ 2 

i 

= d 2 ml { | u x cos 9 cos sin 0' cos 0' + tt y cos 9 sin sin sin 0' — w 2 sin 9 cos 0' | 2 
+ | — -u^ sin sin 9' cos0' + -u^ cos sin sin0'| 2 }. (20) 

Integration over yields 

r-2n 

/ | < m|(e • d)u\l > \ 2 d(f) = -nd 2 ml {\u x \ 2 {l + cos 2 6>) sin 2 9' cos 2 0' 
•/ o 

+ |m j/ | 2 (1 + cos 2 9) sin 2 0' sin 2 0' + 2|m z | 2 sin 2 9 cos 2 0'}, (21) 

which becomes 

/■27T 

/ | <m|(e-d)w|Z> | 2 rf0 = 7r^{|M^| 2 (l + cos 2 ^)sin 2 ^ + 2|M 2 | 2 sin 2 ^cos 2 ^} (22) 

w 

under an assumption of Notice that eq.(22) does not depend on 0' any more. 

When u xy — u z — 1 and 9' = 0, i.e. d is along the z direction, eq.(22) leads to eq.(17), while 
with u xy — u z — 1 and = |, i.e. d is parallel to the x- or y-axis, eq.(22) leads to eq.(18). 
Furthermore, integration of eq.(22) over 9 gives 

/ / | < m\(e • d)u\l > | 2 sin 9d9d<p = nd 2 , \ - \u xy \ 2 sin 2 9' -\ — |m 2 | 2 cos 2 9' \ , (23) 
io io 13 3 J 

which yields §7rd^, when u xy — u z — 1. Actually, eqs.(9) and (10) were obtained by this 
manner. 



2.3 Case of finite space 

In finite space, the electromagnetic field has to be quantized to satisfy the boundary 
condition at the wall. And the energy-momentum tensor is constructed with this quantized 
field which is a solution of the Klein- Fok equation [12]. 

Since the electromagnetic field has a distinctive feature of transverse, the boundary con- 
dition is imposed at the wall such that the tangential components of the electric field vanish 
there. The simplest modification of photon propagator is to discretize the electromagnetic 
field modes in the direction perpendicular to the wall [13]. This photon propagator can be 
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proved[14] to be an approximation of the more reliable photon propagator derived by direct 
insertion of the 5 - function to satisfy the boundary condition[15]. 

All these discussions lead to a result that the energy-momentum tensor gives the dis- 
cretized energy in the vacuum state as 



E= / cfc<0|T 00 |0>=-y> n (24) 
Jo 2^ 

Therefore, the photon emitted spontaneously in finite space must have a discrete energy 

instead of the continuous one. If the hydrogen is placed between two parallel plates of size 

L x L seperated from each other by b(b << L), z-component of the momentum of the 

radiation field which is perpendicular to the plate should be discretized and the integral 

over it has to be replaced by a summation. The result obtained by such a replacement 

is exactly identical with that derived by direct quantization of the electromagnetic field in 

finite space[ll,16,17]. Compare eq.(7a) to eq.(4.7) of ref.ll. Therefore, without repeating 

the field quantization in finite space, we shall start from eq.(7a) to obtain the transition 

rate in finite space. Nevertheless, our theory is rigorously based on the quantization given 

in ref.ll. 

In addition to such a replacement, we must change k z to |n (n = 0, ±1, ±2, • • •) and 

[11] 

-jv u * = \/| exp[-i(k x x + k y y)]cos{k z z) 



exp(— ik ■ r) 



1 



^ exp[— i(k x x + k y y)]i sin(k z z) 



(25) 



where a factor comes from normalization. 

Thus, the transition rate in finite space can be obtained from 

- ~ ?£/«£>i < * ■ <*vw > r^f^. m 

where u x {k z R z ) = u y {k z R z ) = isin(k z R z ) and u z (k z R z ) = cos(k z R z ). R z (0 < R z < |) 
is a distance between the atom and one of these plates. It should be reminded here that 
uj = ck = c(k\ + ky + k 2 z ) l l 2 and the integral over k z is replaced by a summation over n 
because of k z = ^n. Using the relation 
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Ill dk x dk y dk z — = J d(p J d9 sin 9 J dkk = — J d(p J d9 sin 9 J doouo 

— oo 

and replacing the integral over k z by a summation, we find the integral in eq.(26) as 



(27) 



J dVt J duu\ < m\(e ■ d)u(k z R z )\l > 



— / dQ / du\ < m\(e ■ d)u(k z R z )\l > 



,u 3 sin[t(c<Jo — uj)] 

f u 3 sin[t(u — u)] u 3 sin[t(u + u)] 



(ujq - to) 



(UJQ + u) 



7TC 



X 



/oo rco w 
/ dk x dk y ^2 (m|(e- d)u(—R z )\l) 
-oo J -co n=-oo " 

[ fc sin[t{ujQ — ck)] k sin[t(cj + ck)] 1 



(28) 



(c^o — ck) (ljo + ck) J 

where k — [A; 2 + k 2 y + (fn) 2 ] 1 / 2 . By transforming variables, /c x and fcj,, into (ji) and (f j), 
equation (28) becomes 



I = c 3 



7T \ / 7T 



L 



OO /"OO 



— OO J — oo 



(m|(e • eZ)it ( — j^-n) \l) 



x 



Bwa[tu {l-B)] Bsin[tuj (l + B)} 



[l-B) 



[1 + B) 



(29) 



where 



B 



/ C7T \ 2 / C7T \ 2 / C7T 



1/2 



(30) 



Again changing variables as i = {^ ± )x and j = {^ ± )y, we arrive at 



I = zpr r ixdy g 

J -co J -co __ , 



(m|(e • d)u(— -^-n)\l) 



x 



5„ sin[iw (l - 5 n )] £ n sin[tw (l + B„)] 



(1 - B n ) 



(1 + B n ) 



(31) 



where 



-Br) 



2 i 2 



/ C7T 



1/2 



(32) 



In order to evaluate the summation over n, it is convenient to make use of Poisson's sum- 
mation formula on Fourier transforms, 



OO CO ^ 

2 f( n )= / f(u)exp(2msu)du. 

„ _„ ./— 00 



(33) 
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Then, eq.(31) can be expressed as 



7 (±) = £^o jjjdxdydu\(m\{e-d)u(^v^\l) 

x 2 + y 2 + (^ufsm j^o (l ± f^+7+(^f ) } 



exp(27ris-u) 



= - 3 E 



dxdydz 



( m \(e-d)u (—z)\l) 



x 



\/£ 2 + ?/ 2 + -z 2 sin [tuQ (l ± V^ 2 + y' 2 + z 2 ) } ( 2uj bs 



(l ± V^ 2 + V 2 + 



exp 2- 



00 /-oo 



( UqR z 



(m\ (e • d)u r cos 6 |Z) 



c 



s=— 00 

r 3 sin |tw (1 ± r)} f .2uJ bs ^ 
x 7— — ; exp 1 r cos 6> 



(1 ±r) 



Therefore, we obtain 



I = uo, 



00 ;-oo 



(m\(e ■ d)u (^—i-vo>() j /) 



r 3 sin{tcj (1 - r)} r 3 sin{tu> (1 + r)} 



{ (1-r) (1+r) 
3 £ /dft f° dr (m|(e • d)w f^^rcosfl) \l) 



exp %- 



,2u bs 



r cos 9 



s=—oo 
„3 



r 3 sin {tcj (1 - r)} ( .2u bs 
x - ; exp 1 r cos t) 



1X0J. 



where we used 



(1-r) v c 

(m|(e ■ d)w ^ £i!^i_ cos f) j >\ 



00 „ 

e / ^ 



exp 



cos# 



sin[tu;o(l — r)] 



7iS(l — r) 



(34) 



(35) 



(36) 



(1-r) 

for uj t being very large, i.e. c^o = 1-5 x 10 16 s _1 and t ~ 1.6 x 10~ 9 s for the hydrogen 
2P\/2 — > 1 5*1/2 transition. By eq.(22), I can be given as 

. 2 ( ^oRz 



I = n 2 LU 3 d 2 ml £ d9 sin 9 {2 sin 2 cos 0) (l + cos 2 0) sin 2 0' 



+4 cos 2 cos (l-cos 2 #)cos 2 #'} f] exp 



.2cj ^ s 



cos# 



(37) 
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By integrating over 9, eq.(37) is found in the form 

oo 

I = n 2 culd 2 ml ]T {sm 2 ^[4j (6)-2j (ei + 6)-2io(ei-6) 

^'1(6-6 



s=— 00 

4 2 2 

- + 77—7 d i(Zi + 6) + 77 — > v 

+ cos 2 + j^r^hUi + 6) + 7F^TT^(^ - 6)]}, 



(38) 



where j(£) is the spherical Bessel function, £1 = 2uJ0 c Rz and £ 2 = ^f^s. Thus, the transition 
rate is given as 

1 



W = 



8n 3 hc 3 



I = W 1 + W 2 (b), 



(39) 



where 



W x = W { sin 2 9' 



3 3 
1 - 2^0(6) + (£1) 



+ cos 2 9' 



(40) 



W 2 (&) = W sin 2 ^ 



s=l 



3jo(6)-^o(ei + 6)-^o(6-6) 



3 3 3 

+ 9 / f - ; + 6) + 7777 — - 6) 

?2 AU+^J A?l - ?2j 



+ cos 2 0'£ 

s=l 



6 3 3 

t-ji(6) + 77—7 + 6) + 77 — ttJi(6 - 6) 

?2 (,?l+?2j (?l-?2j 



(41) 



W± denotes the s = term, and it is independent of b, while b dependence appears only in 
W2(b). Particularly, when d is polarized along the x or y direction, i.e. 9' — |, W\ leads to 



3 3 
1 - 7>(£i) + 2^Ji(Ci) 



(42) 



If d is along the z direction, i.e. 9' = 0, it becomes 



W lz = Wo 



l + 7-Ji(6) 



(43) 



These results are exactly the same as those given in refs. 11 and 16. 
By defining 

3j'i(¥^-f)) 



Pi = 



:^a-f)) 



(44) 
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3 ,2u; b,l 1 



3jo 
3ji 



2uj b 



2unb 



2 J0 



3ii ( 



s + 



3ji( 



3 
2' 



2u b 



1 Z 
-2 + I 



2(¥0 (*+§-! 



2 ^2o£o&~j 



= E 



s=l 



c 



+ 



)) 



(* + 1 - f) (^) (s-i + !) 

the transition rate can be written as 



(45) 



(46) 
(47) 



W = W [l + sin 2 9'{P 2 + Q 1 (b)} + cos 2 9'{P 1 + Q 2 (6)}] 



(48) 



Here Z = \ - R z (0 < Z < §) is a distance of the atom in the z-direction from the center 
of finite space. 

Since the dipole moment is not polarized in the specific direction in a practical experi- 
ment, we should take average values, < sin 2 6' >=< cos 2 6' >= 1/2, and, then 



W = W {1 + A + A(6)} 



(49) 



with 



A = -(P 1 + P 2 ), 
A(6) = 1[Q 1 (6) + Q 2 (6)]. 



(50) 
(51) 



For b — > oo, the last term in eq.(49) vanishes but the second term, A , can survive as far as 
R z is small, i.e. the atom is located closely to one of the walls. Accordingly, for the case of 
one wall, the modification can be given by A - term. 

§3. NUMERICAL RESULTS 



With various values of b and Z, we can obtain numerical values of the shift of transition 
rate by eq.(49). Our results are shown in Table 1 for the hydrogen 2Pi/ 2 — > 15*1/2 transition 
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with Z = 0. It can be seen that sign of the shift changes in accordance with b value. We 
define 

AW = W -W = W {A + A(b)}, (52) 

values of which are also listed in Table 1. 

For Z = (notice < Z < f ) and b = 1.2/xm, we obtain A = +0.00442, A(6) = +0.0321 
and AW = +22.9MHz, while AW = +38.2MHz for b = 0.7pim. They are in a measureable 
range because the accuracy of the present technique is 4.8MHz[9]. 

The life-time is given by inverse of the transition rate, 

At = t-t (53) 

with r = 1.597ns. Our results are listed in Table 1 and 2. Particularly, Table 2 shows 
dependence of the transition rate on position of the atom, when b = 1.2/im. At Z — 6/10 
which is very close to the center, A + A(6) yields +0.036, i.e. AW = +22.8MHz. This is 
corresponding to At = — 56.1ps. 

§4. CONCLUSION 

In this paper, we have calculated effects of the finite space on the transition rate and the 
life-time of the hydrogen atom. 

We have derived explicitly b dependence of the transition rate, A (b), which was not 
evaluated in the previous calculations [11,17]. The result shows the characteristic oscillation 
of the transition rate in accordance with b values. This property depends very much on 
value of the frequency uo- The magnitude of A(b) varies with value of c^o- 6-dependence is 
around 2 ~ 5% for a separation distance of \xm order. However, it is much enhanced for 
smaller value of u . For example, with u Q = 1.2 x 10 11 Hz, we find that contribution of the 
6-dependence term, A(6), is 33% for b = lcm. 

As is seen in Table 2, variation of the transition rate and life-time are very sensitive to 
position of the atom. Since it is quite hard to set up the atom at an accurate position, 

12 



measurement of the absolute magnitudes of AW and At may be very difficult. However, 
observation of these shifts seems to be definitely possible. 

ACKNOWLEDGMENTS 

This work was supported by the Korean Ministry of Education (Contract no. BSRI-96- 
2425) and the Korean Science and Engineering Foundation. I thank T. Yamazaki and Y. 
Akaishi for their warm hospitality extended at Institute for Nuclear Study, University of 
Tokyo. 



13 



TABLES 

Table 1. Shifts of the transition rate and the life-time for Z — 0. co = 1.5 x 10 16 Hz. 



b(jxm) 


Qi(b) x 10 2 


Q 2 (b) x 10 3 


A x 10 2 


A(6) x 10 2 


AVy(MHz) 


Ar(ps) 


0.5 


-0.1532 


-9.127 


0.038 


-0.5330 


-3.10 


7.94 


0.6 


2.053 


0.9236 


2.423 


1.073 


21.90 


-53.95 


0.7 


9.865 


1.556 


1.081 


5.010 


38.15 


-91.69 


0.8 


-5.800 


1.196 


-1.301 


-2.840 


-25.93 


68.99 


0.9 


-0.5994 


-0.09178 


-1.474 


-0.3043 


-11.14 


28.91 


1.0 


-0.02956 


-1.843 


0.306 


-0.1070 


1.25 


-3.17 


1.1 


1.475 


0.3780 


1.360 


0.7565 


13.26 


-33.10 


1.2 


6.375 


0.5344 


0.442 


3.214 


22.90 


-56.33 


1.3 


-3.090 


0.4390 


-0.923 


-1.523 


-15.32 


40.04 


1.4 


-0.2568 


-0.1551 


-0.858 


-0.1362 


-6.23 


16.04 


1.5 


0.003605 


-0.6429 


0.351 


-0.03034 


2.01 


-4.45 


1.6 


1.281 


0.2144 


0.935 


0.6511 


9.93 


-24.93 


1.7 


4.958 


0.2672 


0.186 


2.493 


16.78 


-41.67 


1.8 


-1.903 


0.2196 


-0.732 


-0.9030 


-10.24 


26.54 


1.9 


-0.1193 


-0.1560 


-0.557 


-0.06743 


-3.91 


10.03 


2.0 


0.02723 


-0.2726 


0.360 


-0.000014 


2.25 


-5.73 
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Table 2. Dependence of the transition rate and the life-time on Z for b = 1.2/im. uj = 



1.5 x 10 16 Hz. 



Zi / 


A \/ 1 n2 
ZAg X 1U 


Za^OJ X 1U 


ZA VV ^IVlllZ J 


A rr- ( -r^n\ 

ZAT^pS J 


U 


n A A o 
U.44Z 


Q 01 /I 

O.Z14 


zz. yu 


-00. OO 


1 /on 
1/zO 


U.841 


O C 1 1 

z.oll 


ZZ.O 1 


ra 0^7 

-00. z/ 


1 / 1 o 

i/iy 


U.4oy 


o.zUz 


OO HQ 

zz.yo 


c^fi /in 
-00. 4U 


i /i c 
l/lo 


-U.Uzoo 


Q /I /1 1 

o.44i 


oi An 
Z1.4U 


-0z. (0 


1/17 

1/17 


-U.004 


o.4o0 


1 7 HQ 

i /.yo 


A A £7 

-44.0 / 


i /i £ 

1/10 


1 n^7£ 
-1.0/0 


3.088 


i o a. i 
Iz.bl 


-ol.ol 


1 / 1 FC 

1/10 


1 /i i n 
-1.41U 


z.o4y 


O.Oo 


1 A 

-14. ob 


1/1/1 
1/14 


1 Q7/1 

-l.oV4 


1 OQO 

i.zoy 


n £qo 
-U.Ooz 


l.ob 


1 /I Q 

l/lo 


-u. I yi 


U.zUoo 


-6.00 


n Qft 

y.oo 


i /in 

1/lz 


O.oOb 


-o.zybo 


n n/^ 
0.0b 


n i £: 
-0.1b 


1/11 
1/11 


l.O ( 


U.4084 


11 A Q 

11.40 


no V£ 
-ZO. ID 


i /in 
1/ 1U 


1.Z00 


z.oVb 


oo on 
zz. OU 


Kft 1 n 
-Ob.lU 


1/9 


-0.618 


3.347 


17.09 


-42.42 


1/8 


-1.478 


0.9838 


-3.10 


7.93 


1/7 


1.528 


0.4776 


12.56 


-31.40 


1/6 


-1.303 


2.974 


10.47 


-26.25 


1/5 


2.088 


1.486 


22.38 


-55.11 
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